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Abstract 

Models of extended gauge mediation, in which large A-terms arise through 
direct messenger-MSSM superpotential couplings, are well-motivated by the dis¬ 
covery of the 125 GeV Higgs. However, since these models are not necessarily 
MFV, the flavor constraints could be stringent. In this paper, we perform the 
first detailed and quantitative study of the flavor violation in these models. To 
facilitate our study, we introduce a new tool called FormFlavor for computing 
precision flavor observables in the general MSSM. We validate FormFlavor and 
our qualitative understanding of the flavor violation in these models by compar¬ 
ing against analytical expressions. Despite being non-MFV, we show that these 
models are protected against the strongest constraints by a special flavor texture, 
which we dub chiral flavor violation (yFV). This results in only mild bounds from 
current experiments, and exciting prospects for experiments in the near future. 
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1 Introduction 


The SUSY flavor problem is a serious challenge for models of weak-scale supersymmetry. 
(For a review and original references, see [1, 2].) Soft SUSY-breaking introduces many 
sources of flavor violation beyond the Standard Model Yukawas. Generic points of the 
MSSM parameter space are ruled out by myriad precision flavor constraints, such as 
neutral meson mixing and b —)■ sy. Clearly, the underlying mechanism that generates 
these soft SUSY-breaking parameters from the hidden sector must be quite special. 
Historically, the lore has been that the mediation of SUSY-breaking must either be 
flavor-blind or obey the minimal flavor-violating (MFV) ansatz. Gauge mediated SUSY 
breaking (GMSB), which is manifestly flavor-blind, is one of the simplest solutions to 
the SUSY flavor problem (see [3] for a review and original references). 

The 2012 discovery of a Standard Model-like Higgs with a mass near 125 GeV [4, 5] 
presents interesting challenges for models of GMSB, especially in the MSSM, where a 
125 GeV Higgs implies either very heavy stops > 10 TeV or large (multi-TeV) stop 
H-terms [6-10]. The heavy stop scenario is more fine-tuned and less interesting from 
both an experimental and theoretical point of view. The large H-term scenario allows 
for light stops, but a mechanism is required to generate these H-terms, which are absent 
at the messenger scale in GMSB. 

Large H-terms can arise if the usual GMSB framework is extended to also include 
direct MSSM-messenger superpotential couplings. Using light (~ 1 TeV) stops and 
large H-terms, these models of extended GMSB (EGMSB) can give rise to the observed 
Higgs mass at fine-tuning levels close to the best achievable within the MSSM [11-28]. 
However, since these MSSM-messenger superpotentials are typically flavorful, they are 
in danger of reintroducing the SUSY flavor problem. Previous works on EGMSB have 
either assumed perfect alignment with the third-generation (to get a large stop H-term), 
or considered additional model building (such as Froggatt-Nielsen mechanisms) to ensure 
this alignment [29-31, 18, 20, 21, 32-34]. In this paper, we will not presume any such 
alignment, and we will instead perform the first comprehensive study of the general 
flavor constraints on EGMSB models for the Higgs mass. 

The precursor to this work was the complete classification of all renormalizable 
EGMSB couplings consistent with perturbative SU{5) unification provided in [20] (see 
also [18]). By turning on one coupling at a time (perfect alignment with the third gener¬ 
ation was assumed) and imposing the Higgs mass constraint, the landscape of EGMSB 
models was surveyed for their phenomenology and fine-tuning. It was shown that the 
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EGMSB models that exhibit the smallest tuning are of the form, 


IV 3 IV 3 K 3 V 3 VV, IV 3 k^Q^HuV, or W 3 (1-1) 

Here <h, $ are messenger helds transforming in appropriate representations of the SM 
gange gronp. The hrst two (second two) models involve a single MSSM held (two MSSM 
helds) and so they were classihed as “type I” (“type 11”) models. 

Notice that the least hne-tnned models in (1.1) are all havorfnl. Type I Higgs models, 
i.e., W 3 kHu^V, while MFV, are more hnely-tnned dne to the “little A/rn\j problem” 
[13]. These models have an irredncible contribntion to + ... where Ai is the 

stop H-term; thns a large Ai is tied to a large and one does not improve tnning in 
these models by increasing the H-terms. 

In this paper, we will stndy the ehects that arise when the conplings of (1.1) are 
no longer reqnired to align perfectly with the third-generation. For simplicity, we will 
specialize to a pair of representative type I Q-class and [/-class models, the models 1.9 
and 1.13 from [20], 


( 1 . 2 ) 


W 3 IV 3 




where the messengers and have the SM gange qnantnm nnmbers of the D and 
L matter helds. We have verihed explicitly that the other type I sqnark models are 
very similar, both qnalitatively and qnantitatively. For reasons we will explain shortly, 
we expect the type H models are also qnalitatively similar in their featnres. We will 
not consider the ehect of tnrning on mnltiple FGMSB conplings, or conplings involv¬ 
ing lepton havor violation. Snch conplings can be forbidden by appropriate choices of 
discrete symmetries. Finally, to focns exclusively on the SUSY havor problem, we will 
not consider GP violation in this work, i.e., all the conplings Ki are taken to be real 
in the mass basis of the standard model particles. The (possibly stringent) constraints 
from GP violating observables, snch as €k and the nentron FDM, will be stndied in a 
forthcoming pnblication [35]. 

The interactions in (1.2) resnlt in havor-violating contribntions to the sqnark mass 
matrices 



(1.3) 


Here each block (clockwise from npper left: LL, LR, RR, RL) is a 3 x 3 matrix. These 
shonld be added to the havor-conserving GMSB contribntions and the snpersymmetric 
contribntions, and together they contribnte to precision havor observables throngh a 
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variety of one-loop diagrams involving the squarks and other superpartners. The calcu¬ 
lation of these one-loop diagrams for general MSSM spectra is a tedious and laborious 
task, ideally suited for a computer program. There are several such programs that are 
publicly available; unfortunately, we found that they were all unsuitable for our purpose. 
These programs either assumed MFV, did not have a sufficiently broad list of flavor ob¬ 
servables, were numerically unstable, or were found to have bugs, likely introduced in 
transcribing formulas by hand from the literature. 

As a result, we found it necessary to develop a new tool called FormFlavor for 
the study of flavor physics. FormFlavor takes general MSSM spectra and computes 
the contributions to the various flavor observables shown in table 1. The novel aspect 
of FormFlavor is that the computation of one-loop Wilson coefficients in the MSSM 
is done completely from scratch, using the general-purpose packages FeynArts [36] and 
FormCalc [37].^ This avoids the problems associated with transcribing formulas from the 
literature, and it facilitates the inclusion of additional flavor observables in an automated 
and modular way. We intend to make FormFlavor publicly available; its usage and 
validation will be described in an upcoming publication [39]. 

Starting from the points of reduced tuning identihed in [20] , we will turn on ki 2 and 
use FormFlavor to investigate the constraints from precision flavor observables. Given 
that they are not MFV, one might expect these constraints to be extremely stringent. 
Generic operator bounds put the scale of flavor violation at 20 PeV or higher from K — K 
and D — D mixing [2]. Even in SUSY, where one benehts from loop factors, etc., the 
bounds are close to 500 TeV (without using GP violation) [49, 50]. However, we will 
hnd that in these EGMSB models, the limits from flavor-violation are extremely mild - 
to the point that ^ 1 , 2,3 can all be the same size and yet the model is not ruled out by 
flavor! 

Much of this paper will be devoted to identifying the reasons for these surprisingly 
mild flavor constraints. Although the overall heavier mass scale required to raise the 
Higgs mass plays a role, the most important reason is the fact that the EGMSB models 
(1.1) only violate flavor through a spurion of either SU{3)q or SU{3)u, but not both. 
Because of this chiral flavor violation (^FV), EGMSB models have a novel flavor texture 
- flavor violation primarily occurs only in either the left-chiral sector of the squark 
mass-squared matrix (for the Q-class models) or the right-chiral sector of the squark 

^While this work was in preparation, a new tool FlavorKit [38] was published with a very similar 
approach to calculating flavor observables from scratch. FlavorKit aims to be even more general, in 
that it can derive the one-loop Wilson coefficients for a general model, not just the MSSM. It would be 
interesting to compare the two codes in detail. 
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Observable 

Experiment 

SM prediction 

ArriK 

(3.484 ± 0.006) X lO”^® GeV 

— 

ArriB^ 

(3.36 ±0.02) X 10-13 QgY 

(3.56 ± 0.60) X 10-13 GeV [40] 

ArriB, 

(1.169 ±0.0014) X 10-11 GeV 

(1.13 ±0.17) X 10-11 GeV [40] 

ArriB 

(6.2+1^) X 10-15 GeV 

— 

Br[K^ —>■ 

(1.7 ± 1.1) X 10-1° 

(7.8 ± 0.8) X 10-11 [41] 

Br{B ^ A, 7 ) 

(3.40 ±0.21) X 10-^ 

(3.15 ±0.23) X 10-^ [42] 

BriB ^ A,7) 

(1.41 ± 0.57) X 10-5 [43, 44] 

(1.54;°;^°) X 10-5 [44] 

Br{Bs -)■ /i+Ai") 

(2.9 ± 0.7) X 10-° [45] 

(3.65 ± 0.23) X 10-° [46] 

Br{Bd -t /U+/i") 

(3.6;j;^) X 10-1° [45] 

(1.06 ± 0.09) X 10-1° [46] 


Table 1: Current experimental values based on PDG and HFAG fits [47, 48] except 
where noted. No reliable theoretical prediction for currently exists. Although 

literature on the subject exists, we do not use a theoretical prediction for Auik (see the 
discussion of ArriK in section 4.2 for more details). 

mass-squared matrix (for the G-class models). Only via communication through the 
SM Yukawas can the other chiral sector feel the flavor violation. As it turns out, the 
most stringent flavor constraints, which come from ArriK and AmD, are vastly reduced 
when the flavor violation is restricted to only the left- or right-chiral sector. (It is also on 
these general grounds that we expect the type II models are similarly unconstrained by 
flavor, although it would be interesting to verify this in detail.) In a forthcoming work 
[51], we will study this new ansatz of X^V in more generality, along the lines of what 
has been done for MFV. X^V represents an interesting intermediate case between full 
flavor anarchy (which is known to be heavily constrained) and MFV (which is known to 
be basically unconstrained). 

In order to validate and interpret our numerical findings, we will compare them 
against analytic expressions for the flavor observables. Historically, the mass insertion 
approximation (MIA) has been utilized to interpret the influence of flavor-violating 
squark masses on precision flavor observables (see e.g., [2] for a review of the MIA 
and original references). However, the utility of the MIA is limited when one or more 
of these mass insertions are 0{1). Since that is precisely the interesting region of 
parameter space for our EGMSB models, the traditional MIA cannot be used here. 
Fortunately, in these EGMSB models there is another handle we can use to obtain an 
analytic understanding. From two powers of our anti-fundamental spurion, Ka, we can 
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construct an adjoint+singlet of SU{3)q or SU{3)u, 

^ab = 


(1.4) 


This matrix K governs all flavor violation from EGMSB. By exploiting the fact that it is 
only rank 1, together with the special properties of the soft masses for the type I models, 
we are able to obtain analytic formulas for the flavor observables that treat the flavor 
violation exactly. Then the only expansion that we do is in v/msusY, which is still an 
excellent approximation in these models, even when the flavor violation is 0{1). With 
this technique, we are able to make precise estimates of the supersymmetric contributions 
to flavor observables, validate our numerical results in detail, and understand their 
qualitative features. 

Our paper is outlined as follows. In section 2, we describe the idea behind chiral 
flavor violation and how it applies to our models (1.1) in restricting the texture of the 
squark mass matrices. We then further restrict our attention to the type I squark models, 
and highlight some of their special features that will be useful in analyzing the flavor 
observables. In section 3, we turn to a detailed study of the flavor constraints on the 
EGMSB models (1.2). We illustrate these constraints in the ^ 1-^2 plane using a series 
of plots of flavor observables computed with FormFlavor. We also provide an analytic 
understanding of the features of these plots using the special properties of rank 1 X^V 
and the type I squark models. Finally in section 4, we conclude with a brief summary of 
our results, and a discussion of the promising future prospects for precision flavor tests 
of these EGMSB models. In the appendices, we detail our parameter deformation and 
various subtleties that arise there, provide a brief description of FormFlavor (postponing 
a more detailed manual and validation for an upcoming work [39]), compile the necessary 
expressions for flavor observables using a uniform notation, and present some other 
formulas used in this work. 


2 Chiral Flavor Violation and EGMSB 

2.1 Chiral Flavor Violation 

In the absence of Yukawa couplings, the Standard Model flavor symmetry group is 

Gf = SU(3)q X SU{3)u X SU{3)d x SU{3)l x SU{3)e (2.1) 

The assumption of MFV is that this flavor symmetry is broken only by the Yukawa 
couplings yu, yd and y£, which transform as ( 3 , 3 ) under SU{3)u x SU{3)q, SU{3)d x 
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SU (3)q, and SU (3 )e x SU{3)l, respectively. As described in the introdnction, the chiral 
flavor violation (xFV) ansatz goes beyond MFV, and postnlates that, in addition to the 
Ynkawas, the flavor symmetry is broken by additional spnrions transforming only nnder 
a single SU{3) of the fnll flavor gronp. 

The ansatz resnlts in a special textnre of flavor violation in the MSSM, one 

which greatly snppresses the constraints from many flavor violating observables. In 
general, the 3x3 soft mass-sqnareds transform as an adjoint+singlet of SU{3)x, 
while the A-terms transform in the same way as the Yukawas. Using spurious from only 
a single SU{3)x, one can obtain flavor violation only in m^. To obtain flavor violation 
in any of the other soft terms, one must involve the Yukawa couplings. Thus these other 
soft masses inherit an MFV-like suppression. 

The cases of interest for this paper are when X = Q or U. (We reserve a more 
general treatment for an upcoming publication [51].) Let us now examine these in more 
detail. To leading order in the Yukawa couplings, the symmetries of Q-class X^V imply 
that 

Smfj = Vu^yi, Smjy = Vd^yl, Au = y^T, Aj = ydV (2.2) 

where S and F are built out of the SU{3)q spurious. Substituting these into (1.3), 
we have under the third-generation dominant approximation (whereby yu and yd are 
nonzero only in the 33 component): 



/ 




0 

0 

mtF*i 

\ 




5ml 


0 

0 

miF*2 


6Ml^ 





0 

0 

mtF*3 



0 

0 

0 

0 

0 

0 




0 

0 

0 

0 

0 

0 



\ 

31 

mj:32 

mt^33 

0 

0 

y^^33 

/ 


/ 




0 

0 

m^Tgi 





^rriQ 


0 

0 

rUbT ^2 


SMl ^ 

a 





0 

0 

m^Fgg 



0 

0 

0 

0 

0 

0 




0 

0 

0 

0 

0 

0 



1 

mbTsi 

mbT32 

rribTss 

0 

0 

y^b^33 

J 


We see that within the third-generation dominant approximation, no flavor violation 
appears in the RR block for either the up or down squark mass matrices. Also, the 
only flavor violation in the LR block involves the 3rd generation, and is v/msusY sup¬ 
pressed. These features greatly reduce the sensitivity of Q-class X^V to precision flavor 










constraints. As we will explain in more detail in section 3.1.1, what typically provide 
the most stringent flavor bounds on new physics (i.e., ArriK and Amo) involve the 1st 
and 2nd generations and involve simultaneous violation of flavor in both the left and 
right chiral sectors. 

Meanwhile the [/-class models are even more insulated from constraints. Here, 
on symmetry grounds, and again to leading order in the Yukawas, we must have 

= vi^Vu, As = Tyu, = Aj = 0 (2.4) 


where now S and F are built out of the SU{3)u spurious. Again substituting these into 
(1.3) in the third-generation dominant approximation, we And: 



/ 

0 

0 

0 

0 

0 

0 

\ 


/ 







0 

0 

0 

0 

0 

0 



0 

0 

0 


SMl ^ 


0 

0 

y^^33 

mtF*3 

to * 

CO 





0 

0 

0 

0 


0 

0 

13 






0 

0 

y?^33 




0 

0 

mt^23 


5mlf 




V 



0 

0 


\ 

0 

0 












Thus there is no down-type flavor violation at all in the third-generation dominant 
approximation! As the majority of sensitive flavor probes involve the down sector, this 
eliminates almost all contributions to flavor observables. The only exception is Amo, 
which again has reduced sensitivity because there is no simultaneous left-right flavor 
violation in the lst/2nd generations. (Down-type flavor observables can involve the up- 
squark RR block through chargino loops, but these must be mostly Higgsino-like, so 
they will be suppressed by Yukawa couplings.) 

Finally, we should comment on the role of the RG. The RGBs from the messenger 
scale to the weak scale add more terms, but the symmetry-based arguments given above 
- which were truncated at leading order in the Yukawas - clearly continue to hold with 
the inclusion of higher orders in the Yukawas. In particular, terms that are zero in the 
third-generation dominant limit remain zero under RG evolution. For this reason, we 
do not need to concern ourselves with the details of RG running to gain a qualitative 
understanding of the EGMSB flavor-violating contributions to our flavor observables. 

2.2 aFV in Type I Squark Models 

The EGMSB models (1.1) are clearly examples of Q-class or [/-class XFV. In this sub¬ 
section, we will discuss some further features that are specific to the type I squark models 


9 










that will help us analyze the Wilson coefficients for flavor observables in the next section. 

We begin by quoting the explicit formulas for the soft masses for general type I 
squark models (see appendix C). For concreteness, we focus on Q-class models, where 
we have (in addition to the GMSB contributions). 


^^2 _ dQ 

~ 2567r4 


2 ^2 f A\ A? 


(dfj, + dQ)K — 2Crgj. — 




2 _ dqdjj 

“ 2567r4 


5rn?H = 
= 




2567r4 

dqA 


yuKy\A^ 

Ti[yuKyi]A^ 


Ibvr^ 


yuK. 


( 2 . 6 ) 


Here we are neglecting the (numerically irrelevant) down-Yukawa contributions; K is 
the rank one matrix of couplings dehned in (1.4); d^^ = 2 is a multiplicity factor; and 
h{x) is an O (1) loop function (see (C.8) for the exact form). For our Q-class model of 
study (1.2), we have dq = N = 6, d^p = 5, where N is the number of messengers.^ 

We can see the texture of (2.3) quite clearly in the formulas for Smfj and Au in 
(2.6). There are also additional features of these explicit formulas that go beyond the 
XFV ansatz. We notice that Srriq, together with the quantities F and S introduced in 
(2.2), are all proportional to K: 


5ml = /3K, F = yiF, E = aK 


(2.7) 


The simplicity of these relations is partly due to the rank 1 nature of the flavor viola¬ 
tion. But in principle, on symmetry grounds alone, there could have been additional 
contributions to (2.7) proportional to the identity matrix and powers of the Yukawa 
couplings. These are absent due to the specihc form of the type I squark couplings. 
This is the main reason we have focused on the type I models in this paper. The forms 
of the soft masses are more complicated in the type II models, and while we expect them 
to be similarly protected by their X^V flavor texture, understanding them at the level 
of analytical detail that we apply to the type I models is more difficult. 

Using the relations (2.7), we now discuss the diagonalization of the squark mass ma¬ 
trices. Because the LR blocks are suppressed by v/msusY, the squark mass eigenvalues 
are given by those of the LL and RR blocks to a very good approximation. The RR 

^For the similar {/-class model formulas, simply change U -O- yq —)• j/|, and K —)• K'^. 

The analogous multiplicity factor is = 1, and for our {/-class model, we have du = 2Y = 6, d(l> = 4- 
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blocks are already diagonal in the third-generation dominant approximation according 
to (2.3). The down-squark RR masses are just their minimal GMSB values ttIq, while 
the up-squark RR masses are given by mg, ml and mg -|- crylnl = m\j^. In practice the 
third eigenvalue is only a little offset from mg throughout the parameter space of our 
model, and this difference can be neglected. 

The situation for the LL block is only a little more complicated. Since ^mg is just 
proportional to K, the LL block of the squark mass matrices is diagonalized by the 
unitary transformation 



( ki 

Oi 

bi\ 

Ull = 

k2 

02 

b2 


yks 

03 

bsj 


where ka = Kq/k and a and b are any two orthonormal basis vectors for the orthogonal 
subspace /?■*“. The eigenvalues are mg, mg, and 

m| = mg -1- / 3 k ^ (2.9) 

(Here we are ignoring the small differences between the minimal GMSB contributions 
to the Q, U, D soft masses coming from SU{2) x G(l).) In our EGMSB models, m| 
corresponds roughly to the mass of the lightest squark, and moreover 

m| -C mg (2.10) 

Thus the Wilson coefficients will generally be dominated by this lightest squark running 
in the loop, and m| will play a central role in controlling the size of the flavor-violating 
effects in these models. In fact, because of (2.10), diagrams where the right-handed 
squarks propagate tend to be suppressed, and we will see that it suffices to focus on the 
LL flavor violation exclusively. 

3 Flavor Observables and Constraints on EGMSB 

In this section, we investigate the flavor constraints on our EGMSB models (1.2) in 
detail. For reference, the spectra corresponding to = ^2 = 0 are shown in fig. 1. 
These are essentially the same points that were identified in [20] as being the least fine- 
tuned EGMSB models with mh = 125 GeV. Starting from these flavor-aligned points, 
we will perform a numerical scan in the flavor-violating parameter space {ki/k^, K 2 /K 3 ) 
of the models, using FormFlavor to compute the flavor observables and compare against 
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their experimental values. Going through each observable X from table 1 in turn, we 
will exhibit plots of 


[■y]T„ - [Vexp 

TH+EXP 


Here denotes the SM prediction together with the EGMSB contribution, and 

is the theoretical and experimental errors added in quadrature. A contour 
of 2, which is roughly the 95% GL exclusion limit, will be taken to indicate the point at 
which the EGMSB model is excluded by the given flavor observable.^ 

As we deform away from ki = K 2 = 0, there are numerous subtleties that must be 
taken into account regarding how the other parameters of the model are varied. These 
subtleties and specihcs of the procedure are described in appendix A. In short, ki and K 2 
are introduced in such a way that the superpartner mass eigenvalues and “net” A-terms 
are essentially held hxed. 

In order to validate the numerical results from FormFlavor, we will compare them 
against analytical formulas for the flavor observables. This will also shed further qual¬ 
itative insights on the role of weakening the flavor constraints. As discussed 

in the Introduction, the usual mass insertion approximation fails due to O {1) entries. 
Instead, we will use flavor symmetries and the special features of the type I models dis¬ 
cussed above to characterize their exact k dependence, to leading order in an expansion 
in v/ms- In general, the k dependence arises through the lightest squark mass (2.9) 
and the unitary matrix Ull given in (2.8). (The LR and RL blocks in (2.3) also depend 
on K through (2.7), but as discussed in the previous section, this dependence can be 
generally be neglected due to the heaviness of the right-handed squarks and the extra 
v/ms suppression.) The dependence on Ull is constrained by the flavor symmetries, 
and as explained in the previous paragraph, the squark mass eigenvalues are mostly held 
hxed in our parameter space. Thus it will be possible to fully characterize the features 
of the FormFlavor plots in terms of very simple functions of k . 

^Two exceptions to this are ArriK and Aitid where the theory errors are uncontrolled. In these 
cases, we will plot instead: 


I [ATOx]EQ]y[3g I 

[Amx + 2(T(Amjf)]g^p 


(3.2) 


Thus a contour of 2 does not represent a 95% confidence level exclusion for these observables. However, 
as it would require a substantial tuning for the standard model and the new physics contributions to 
cancel against one another, values larger than 1 are suspect. 
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Figure 1: Top: Mass spectrum of the model at ki = K2 = 0. Bottom: Mass 

spectrum of the KiU ^2 model at ki = ^2 = 0. 
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3.1 Q-class Models 

3.1.1 Meson Mixing 

We begin with the AF = 2 meson mixing observables. As we discuss in appendix D.l, 
for X = K and D, 

Amx = 2 Re (X\Heff\X) (3.3) 

to a good approximation, while for X = Bd and Bq, 

l^mx=2\(X\H,„\X)\ (3.4) 

to a good approximation. Here He// is the AF = 2 effective Hamiltonian; its local 
short-distance part (relevant for the MSSM contributions) is built out of the four-fermi 
operators, 

(CO'").! = (iiaYPMMq.l^PNqi) (3.5) 

Here M, N = L,R label the chirality of the incoming quarks; and ab = 12,12,13, 23 for 
ArriK, Amo, Attib^ and Auib^, respectively, while q is up-type for Auid and down-type 
for the rest. 

The full result of FormFlavor is shown in hgure 2. In this subsection, we will endeavor 
to understand its features analytically using the X^V ansatz and the special features of 
type I EGMSB identihed in section 2.2. 

The MSSM contributions to the AF = 2 observables are due to box diagrams in¬ 
volving the squarks and the gauginos. Here the great simplihcation of X^V is that any 
operator with an R index must transform non-trivially under SU{3)b> or SU{3)u] thus it 
is suppressed by X^V in our Q-class models. Furthermore, the and operators 
all involve an SU{2) ^-breaking chirality flip, so they are dropped in the n = 0 approxima¬ 
tion. Therefore, the only unsuppressed Wilson coefficient is Cy^. This also happens to 
be the only contribution to the one-loop SM Wilson coefficients, which proceeds through 
W exchange. 

Since Cy^ transforms in the square of the adjoint-l-singlet representation of SU{3)q, 
with just left-handed squarks running in the loop, the only way it can depend on k is, 

(G“)., = ^/i(m|/m2), (3.6) 

TUs 
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Figure 2: Plots of meson mixing observables the Q-class model, KiQi^D^L- 
Amo, Amo^, Attib^ are presented in the upper left, upper right, lower left and lower 
right, respectively. Both Am^ and Amo are presented as the FormFlavor output over 
the experimental value, while Am^^ and Am^^ are the difference between FormFlavor 
and the experimental value in units of the uncertainty. 

where /i is a dimensionless loop function depending on the LL squark mass eigenvalues.^ 
Dimensional analysis hxes the dependence on the masses, and the rest of the dependence 
must be from the unitary matrix (2.8). Under this simplihcation, the meson mixing 

■^All loop functions here and below will implicitly depend on the masses of the other superpartners 
running in the loop, e.g., the gaugino masses. 
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contributions are of the form (see appendix D.l for more details and specific values of 
the parameters), 

\(X\H,,,\X)] 

EGMSB ~ ^'rnxfx Bv!x{^V^)ab, (3.7) 

where fx is the decay constant for the meson; and By^ is an (9 (1) hadronic parameter. 
This simple formula suffices to accurately describe the AF = 2 flavor observables. 

Focusing on gluino boxes for simplicity,® we hnd from explicit computation that the 
leading Wilson coefficient is given by 


{Cy^)ab — kl 


a. 


-ft 

J a 




[Xg,Xq) 


(3.8) 


This clearly agrees with the general form (3.6). Here, Xq = mg/m^ and Xg = m'^/mQ, 
and Xq) ~ 0.05 is a loop function that we dehne in appendix E. As discussed 

earlier, since the relevant masses are roughly held hxed in our deformation, the loop 
function does not change across the parameter space. The QCD RG running from the 
SUSY scale to the meson scale is fairly mild for Cy^, yielding only a 20 — 30% suppression 
in the size of the Wilson coefficient [52]. 

The hadronic factors fj^By\ are fairly similar across the mesons. We can dehne. 


Hx = 


f. 


2 txLL 

x-‘->v,x . 

f2 nLL 5 

J kx>v,k 


Hk = F H 


D 


H 




2 , H 


Bs 


(3.9) 


So we expect deviations of the form. 


[(X|77e//|X)]EGMSB ~ -kl klHxmx (lO-^®) , (3.10) 

Noting that kf^kl is at most |, and comparing against table 1, we see that Amx and 
Amo should be most sensitive to EGMSB, while ArriB^ should be barely sensitive, and 
ArriBs completely insensitive.® 

These sensitivities are observed in the plots shown in hgure 2. Moreover, the k^ k^ 
dependence is transparent in these plots. AB^ peaks at (ki/ks, ^ 2 /^ 3 ) ~ (il^O); while 

®Chargino and gluino-neutralino boxes are also typically comparable. However, they do not affect 
the qualitative discussion here. In fact, because the charginos enter with the opposite sign of the gluino 
and gluino+neutralino contributions, the gluinos alone provide a better estimate quantitatively than 
might be expected. 

®As discussed in the Introduction, we are not considering CP violation in this work, in particular 
Efc- Although it depends on the precise value of the CP violating phase, the expectation is that this 
will place a meaningful, tighter constraint on the parameter space when the phase is large. This and 
other CPV observables will be studied in an upcoming paper [35] . 
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ABg peaks at ~ (0, ±1). For ArriK and Am^i, we can see that moving from ki = K 2 = 
to the corners, constraints rise by the expected factor of ~ 16/9/ 

These results from meson mixing may seem to conflict with the SUSY flavor problem, 
which suggests that, with O (1) flavor-violation, the SUSY scale needs to enter above 
~ 500 TeV due to constraints from Attik and Amo [49, 50]. However, these constraints 
are driven by the Wilson coefficient, while Y^V only generates the Cy^ operator. 
In the MSSM, the contribution of the latter to AmK,D is suppressed by 10“^-10-3 
relative to the former. This is due to three separate effects that all work in the same 
direction. First, the hadronic matrix elements differ between these two operators. 


{K\Slr\K) 

{K\Vll\K) 


A~B^ 


~ 35, 


{D\Slr\D) _ 3Bl^ 

(BIVilIB) 4S“ 


(3.11) 


Next, the SUSY contributions to the Wilson coefficients are also quite different in size. 
From the MIA with O {!) mass-insertions, it is easy to see that Cg^/Cy^ ~ 30 [49]. 
Lastly, the QCD running from the SUSY scale to ~ 2 GeV suppresses Cy^ by 20-30%, 
but enhances Cg^ by a factor of ~ 3 [52] . All of these factors conspire to drop the scale 
of sensitivity to Amx and Amo to the TeV scale in models. 


3.1.2 

Unlike meson mixing in the previous subsection, the —)■ observable (along with 

all other AF = 1 observables) enters as a matrix element squared, and thus interfer¬ 
ence with the standard model contribution can be important. The expression for this 
branching ratio is (see appendix D.2), 

BR(iF^ ^ TT^uu) = c+v^\C^^gj^ + (3.12) 

where c+v^ = 4.9 x 10^ GeV^ C^jg^ = (-1-21 + 0.39i) x 10-3° GeV-^ and all EGMSB 
effects are contained in Cy^ and Cy^. These are the Wilson coefficients from the AF = 1 
effective Hamiltonian built out of the four-fermi operators 

(Of "(ai = (3.13) 

(Of").! = 

^Note that Amo is rotated by the Cabbibo angle G relative to the other observables. This is a 
consequence of the the fact that, as discussed in in Appendix A, our ki, K 2 , and K 3 directions are 
chosen to align with the low energy down, strange and bottom quark. For Amo, the dependence on Ki 
proceeds through the LL sector of the up-squark mass matrix, where there is an additional rotation by 
VcKM- 
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Figure 3: Contours of the difference in BR(iC’'' —)■ n^uV) between experiment and the 
SM + EGMSB predictions, given in units of the net uncertainty (added in quadrature). 
A contour of 2 corresponds roughly to a 95% exclusion. 

Here ah = 12 for —>■ In the MSSM, these Wilson coefficients arise through 

one-loop box and Z-penguin diagrams. 

The full FormFlavor result is shown in hgure 3. The general trend of —)■ 

can again be understood through use of the features of rank 1 discussed in section 
2.2. First of all, as for the AF = 2 observables, must be zero in the third-generation 
dominant approximation - since it transforms in the adjoint-l-singlet of SU{?))d^ using 
the available spurions one can only obtain something that is nonzero in the 33 compo¬ 
nent. Thus, we can focus on Cy^f^. This transforms in the adjoint-Fsinglet of SU{3)q. 
Using flavor violation in the LL block only, the form of is constrained by the 

symmetries to be: 



(3.14) 


where ... contains higher orders in v/msusY and other irrelevant terms, and /2 is a 
dimensionless function of superpartner mass ratios. We note that (3.14) can come from 
box diagrams or Z-penguin diagrams, but in the latter case, the l/ml from the Z 
propagator must be canceled out by two insertions of wino-Higgsino mixing. Insertions 


of LR mixing from the squark mass matrix would also cancel out the l/ml, but the 


heavy right-handed squark masses suppress these contributions enough that they may 
be ignored. 
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Now we will compare against explicit computations of the Wilson coefficients. We 
confirm that the Cy^ coefficients are all negligible. For Cy^, we hnd that the chargino 
diagrams dominate, and they are given by 

(3.15) 

Again, ~ 1.5 is a loop function of the sparticle mass-ratios = 

X 

mym'g, X 2 = and x^ = jji^/rn?g, dehned in appendix E, and varies only mildly 

across the parameter space. This is fully consistent with the general form (3.14). This 
also explains why the charginos dominate over the gluinos, since the gluinos can only give 
rise to Z-penguin diagrams that are nonzero by virtue of down sector LR squark-mixing 
insertions. 

Substituting in numerically for the loop function, irts and 02 , we hnd 

~ (8 X 10"^^ GeV^) kik2 (3.16) 

This simple function of k, when added to the SM contribution and substituted into 
(3.12), reproduces well the features of hg. 3. It grows in magnitude fastest along the 
lines Ki = ±^ 2 , asymptoting to the values ±4 x 10“^^ GeV“^. In the corners of the 
parameter space we obtain a deviation from the SM prediction of 

ABR(A:+ ^ 71+uu) = c+v^ (^2Re [C^^smCv^*] + ^ ^5 x 10-^\ (3.17) 

This is smaller than the current experimental-|-theoretical uncertainty shown in table 1. 
However, because the SM prediction for the BR is a little lower than the experimentally 
observed value, moving along the ki = +K 2 line slightly exacerbates the difference, while 
moving along the Ki = —K 2 line slightly lessens it. 

3.1.3 b ^ S'j and 5 —)■ dy 

The expression for the b ^ sj branching ratio (see appendix D.3) is, 

BR(6 ^ S7) = + Ci|^ + |C«|^), (3.18) 

where = 1.97 x lO^^ QeV^ and C^sm ~ 1-3 x 10"® GeV^ The new physics is 
contained in and C^; these come from the effective Hamiltonian built out of the 
dimension 5 operators 

{0^)ab = e{q,a^’'PMqb)F^, (3.19) 
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Figure 4: Sensitivity to 6 —>■ sy (left) and 6 —)■ dy (right), same conventions as fig. 3. 


where ab = 32 for b —)■ sy. These operators require SU{2) breaking, so their Wilson 
coefficients are zero in the n —)■ 0 limit. Naively, this would mean that the result is 
negligible in our zeroth-order expansion in v/msusY- However, the SM Wilson coefficient 
is suppressed by rrib/v, so here we consider one higher order in the v/msusY expansion 
in order to capture a numerically relevant result. 

Since (C^) transforms in the (3,3) of SU{3)£) x SU{3)q {SU{3)q x SU{3)d), it 
must be proportional to rud acting on the left (right). The latter is zero in the third- 
generation dominant approximation, so it suffices to focus on With just left-handed 
squarks propagating in the loop, the Wilson coefficient must be given by: 


_ {mdK inn ( 3 ) 2.2 1 ^ ^ ^ rubtan/^ 

-zry/e - —zy —h 




m% 


(3.20) 

og Trig 

Here we assumed that rrid is accompanied by a tan/? enhancement, otherwise the entire 
effect is numerically negligible. 

Comparing with explicit calculation, we hnd again that charginos give the dominant 
contribution to the Wilson coefficient, through the quark-squark-Higgsino vertex. (Con¬ 
tributions of the form (3.20) can also arise through gluino and neutralino loops, but here 
the factor of mb tan (3 arises through the LR block of the down-squark mass matrix, so 
the diagrams are suppressed by heavy right-handed squarks propagating in the loop.) 
The dominant chargino diagram gives. 


= K K ^^*^2 pb^s/d'y,pens/ N 

^ 2^3 9 ./v± jJ.) ^2)) 


mi 28871 ^ 


(3.21) 
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where = fi^/rrig, X 2 = M^/rrig, and ~ 0.5 is a loop function 

defined in appendix E. This is fully in agreement with the general result (3.20). QCD 
running from the SUSY scale to the b pole [53] induces a mild 20% suppression to (3.21), 
and once again, we find this simple result is enough to account for the features of the 
FormFlavor plot in figure 4. 

Substituting in numerically for our parameter space (we take tan/? = 10 as in [20]), 
we obtain: 

Ui ~ (3 X 10“® GeV^) kska (3.22) 

Clearly its effects are largest when ki = 0, and along the ^2 axis it has a maximum 
(minimum) at ^2/% = 1 (-1), while asymptoting back to zero as k,2/k,3 —)■ ±00. Substi¬ 
tuting this into (3.18) together with the SM Wilson coefficient, we find at these extrema 
a deviation from the SM prediction of 

[ABR(6 ^ S7)1...±.3 ~ c-r”" (2 [CIsmC\’] ) ~ ± 7 x 10-^ (3.23) 

Since the uncertainty on the measurement (combining the theoretical and experimental 
in quadrature) is 3.5 x 10“®, regions of exclusion are to be expected for h —)■ sy. 

In practice, FormFlavor does hnd a region of exclusion in the bottom half of the 
parameter space, shown in hgure 4. This exclusion is in part due to the current O (la) 
excess in the experimental measurement relative to the theoretical prediction,® which 
allows the ~ 2a change to constrain a large region of negative K 2 in the plot. Uncertainty 
on this observable is comparable in size between theory and experiment, so improvements 
on either side could make this observable more constraining and in a direction that other 
observables currently have no sensitivity. 

The observable 6 —)■ dy could also potentially place constraints. The EGMSB con¬ 
tribution to the Wilson coefficient is the same as in (3.21)-(3.22), but with K2 O ki. 
However, the SM contribution is of a different size, ~ — (2.4 -|- 1.1?) x 10“® GeV“^ 

Importantly, this is not much larger than the new physics contribution, so interference 
is very important. This yields an approximate deviation of, 

lABR(t ^ £(7)],.=., (-2 |Re | + |C‘f) ~ -0.8 x lO'^ 

^ ^ (3.24) 

[ABR(6 ^ <i7)l„._.. (2 |Re | + [Cif) ~ 1.6 x 10-^ 

^Recently, an improved theoretical prediction of 6 —> sy was released [54]. This work predicts that 
BR(6 —)■ sy) = (3.36 ± 0.23) x 10“^, which is significantly more in line with the experimental value. 
Unfortunately, insufficient details are provided in that work for us to modify FormFlavor to account 
for this prediction. 
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Figure 5: Sensitivity to Bg —)■ /i/r and —)■ /x/i, same conventions as fig. 3. Neither is 

constraining over the parameter space. 

while the net uncertainty is 6.3 x 10“®. Again, a two sigma exclusion could be possible, 
but it would only be expected near ^ = — 1 due to the interference acting constructively. 
This is borne out by figure 4, which shows a maximum deviation of ~ 2.3 times the net 
uncertainty. Because the measurement of 6 —)■ dy is fairly recent, a future improvement 
that pushes the experimental uncertainty to the 10-15% level would allow for 6 —)■ dy to 
place much tighter constraints on the region of constructive interference, i.e., Ki < 0. 

3.1.4 Bq —>■ /r ■*■/!“ 

As shown in £g. 5, neither Bg —)■ nor B^ —)■ is at all close to constraining our 

EGMSB models. We include a brief discussion of these observables just for completeness 
sake. 

The branching ratio for Bq —)■ is [55, 56] (see appendix D.4 for more detailed 

formulas): 

BR{Bq ^ p+/i-) = Xg I 1^1 - ^ j (3.25) 
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Figure 6: Wave function correction diagrams that provide the leading, tan^/9 contribu¬ 
tions to Bg —)■ Loop can be a gluino, chargino, or neutralino. 


where Xg = 5.36 x 10^ and Xd = 3.97 x 10^, and 

Q 

7T7 

171 ( 9 ) _ c/^LL I /^LR /^RR /^RL\ 

“ + 

Til 

^ mfe + m/ * * * * 

- C^^), 


Here the Cy^ coefficients are those for the AF = 1 effective Hamiltonian introduced in 
(3.13). The three-loop Standard Model contribution [57] is F^'^gj^ = (1.5 — 0.6i) x 10“® 
and F^^gj^ = (—7.9 — O.lf) x 10“®. 

In the MSSM, the dominant contributions to Bq —>■ are the tan^ (3 enhanced 

wave function correction diagrams with a heavy higgs propagator, see hgure 6. Here 
gluinos, charginos and neutralinos can run in the loop, with both CP even and odd 
higgs states along the penguin line. For simplicity, we will quote the result only for 
gluinos; the answer for the others is very similar. As our higgs states are heavy, we use 
the relation m\ to simplify expressions. 

These diagrams contribute to , and the discussion is similar to that of 6 —)■ gy. 
Cg^ transforms in the (3, 3) of SU{3)d x SU{3)q while Cg^ transforms in the (3,3) of 
SU{3)q X SU{3)d- Thus as before, the latter is zero thanks to X^V and we can focus 
on the former. Here we can afford to work at n = 0 since we are not concerned with 
Z-penguins. With only left-handed squarks propagating in the loops, by symmetries the 
answer must be of the form: 


yyLM _ 
'-^5,3a — 


{.^Vdjsa p _ ^ ^ Vb p 

-Jb — -y/s 




(3.27) 

og Trig 

where a = 2 for —)■ fifi and a = 1 for —)■ /i/r. Now the Yukawa coupling needed on 
symmetry grounds arises from the Higgs-quark-quark coupling. Explicit computation 
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gives 


^LR 




mb tan (3 m^ tan^ (3 Aa 2 a 

^ o 


mw 


mw 3m\ ms 


f 


H ' r~ 


^g) 


(3.28) 


where Xg = m‘^/ml and Xq = m|/mQ and the loop function (see appendix E) is 
Xg) ~ 0.1. Since m^ ~ 'msusY: this leading contribution is of the 
form (3.27) as expected. 

From the above expressions, we can translate to the phenomenologically useful pa¬ 
rameters 


Fjf) = ~ (2 X 10-^°)/t3K2 

^ mb + m^ ^ ^ ^ 

p(d) ^ ^lr _ ^2 X 10-^°)/t3/ti. 

mb + ma 


(3.29) 


These translate into maximum deviations from the SM branching fractions by 


[ABR(i?5 —)■ /i^/r’ 
[ABR(i?d —!• 


MAX 


MAX 






2 Re 
2 Re 


p(®)* p(^) 

^ A,SM^ P 

T?{d) 

^ A,SM^ A 


~ 9 X 10"^^ 

~ 1.2 X 10"^^ 


(3.30) 


These deviations are an order of magnitude smaller than the uncertainties on their 
respective measurements, and therefore neither Bg —)■ nor —?• fi^/ jT place any 

meaningful constraints on our models.® 


3.2 U-class Models 

As discussed in section 2.1, all 17-class X^V models receive very few constraints from 
flavor observables. First, according to the texture of [/-class models (2.5), all flavor- 
violation is restricted to the up-squark sector only. As most potentially constraining 
flavor observables have external down-type quarks, chargino diagrams are required for 
sensitivity to flavor-violation. However, as is further shown in (2.5), there is no flavor- 
violation in the LL block; thus pure wino diagrams cannot contribute. Higgsino diagrams 
introduce Yukawa couplings which suppress contributions to down sector observables 
enough that none of these would be remotely constraining in the foreseeable future. We 
have verified all of these general results in the context of our type I EGMSB models. 

There can be constraints from the //-meson system. Contributions to Am^ are as in 
the Q-class models, only now with the Cy^ contribution dominating, so we expect very 

®In order for Bg -A to place any meaningful constraint, one would have to take tan/3 much 

larger than the tan/3 ~ 10 that we have assumed in this paper, e.g., tan/3 > 30. 
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Figure 7: Plot of Amn meson mixing observable for a fZ-class model. 


similar constraints from Am^). This is indeed shown in hgure 7. Again, although the 
contour of 2 does not represent a 95% conhdence level exclusion, values larger than 1 
necessitate a cancellation between the standard model and the new physics contributions. 
It should be noted that a viable possibility is that the standard model contribution is 
in fact much smaller than the observed value, and the contour of one is actually where 
EGMSB entirely accounts for the ArriD measurement. As with the Q-class models, the 
sensitivity vanishes as either ki —)■ 0 or ^2 —>■ 0, and increases most rapidly along the 
diagonals. 

4 Conclusions 

4.1 Summary and Discussion 

In this work, we performed a detailed investigation into the precision flavor constraints 
on extended GMSB models. These models, where the usual GMSB contributions are 
augmented by direct matter-messenger couplings, are well motivated in light of the 
recent discovery of a Higgs boson near 125 GeV. However, since these models are not 
necessarily MFV, in their full, three-family generalizations they could be potentially 
dangerous from the point of view of flavor. 

Our work required a computer program that could turn general MSSM spectra into 
precision flavor observables. We found that existing programs had various limitations 
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- they either assumed MFV, were numerically unstable, or had incorrectly transcribed 
formulas from the literature. This motivated us to develop FormFlavor, a comprehensive 
package that computes flavor observables ah initio starting from the Feynman rules, and 
uses a modular framework that enables us to add new observables in a uniform and 
straightforward way. 

Using FormFlavor, we studied the flavor constraints on the three-family generaliza¬ 
tions of the EGMSB models of [20]. The results we encountered from this systematic 
study were interesting and unexpected. We found that despite the introduction of O (1) 
flavor-violating couplings, there are currently very few constraints on EGMSB models. 
To validate the numerical results of FormFlavor, we compared them in detail with ana¬ 
lytic formulas for the Wilson coefficients derived using a combination of flavor symmetry 
arguments and direct calculation, and we found excellent agreement. 

The mild flavor constraints in these models are illustrated in the summary plot 
of hgure 8. [/-class models only receive constraints from Am^. Q-class models are 
constrained by Atuk and Atud in the corners of the plot, while the radiative b ^ s'j 
and b ^ d'j each exclude a single bubble near k ,2 = — 1^3 and ki = —^3, respectively. We 
note that the excluded region from 6 —)■ sy is largely due to the current ~ Icr discrepancy 
between the theoretical prediction and the measurement, which new theoretical work 
suggests will disappear [54]. 

The results from these models may seem at odds with the SUSY flavor problem. We 
have argued that the mildness of the flavor violation in these models originates from the 
fact that they obey the “chiral flavor violation” ansatz, whereby flavor is violated only 
by the Yukawas and spurious of a single SU{3) of the full SU{3)^ SM flavor symmetry. 
We showed that X^V prevents many of the most problematic contributions to flavor 
observables from arising in the MSSM, and allows for O (1) flavor-violation in EGMSB 
models. 

4.2 Future Constraints 

Although constraints are currently very mild, there is immense potential to further probe 
these models in the near future. 

• For AmK, short-distance predictions exist with moderate uncertainty, i.e., (3.1 ± 
1.2) X 10“^® GeV [58], but the long-distance contributions are currently unknown. 
However, an accurate, full calculation on the lattice, including both long- and 
short-distance contributions, may be coming in the near future, as promising pre¬ 
liminary work on the subject shows [59, 60]. 
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Figure 8: Flavor sensitivity for a Q-class (left) and ?7-class (right) model. Shaded 
regions of this EGMSB parameter space are excluded to 95% conhdence level by current 
measurements of flavor observables as described in the text. Dashed lines indicate the 
locations where the EGMSB contribution to ArriK and Amn are equal to the measured 
value (as these currently cannot provide a genuine exclusion due to poorly controlled 
theoretical uncertainties). In both cases, observables not shown are not constraining. 

• Expected incremental improvements to the theoretical uncertainty of the standard 

model prediction could make AB^ into a constraining observable soon. Estimates 
suggest that both the bag parameter, and the relevant GKM elements could 

be calculated on the lattice to signihcantly improved levels by 2018 [61]. 

• Belle II [62] is expected to make signihcantly improved measurements to both 
6 —)■ S 7 and 6 —)■ dy [63], which would allow for both of these observables to 
constrain more of the parameter space. 

• NA62 [64] at GERN is projected to be able to measure BR(iF=*" —)■ 71^ uu) to a 
precision of about 10% [65] of the SM value. This is an improvement of more 
than an order of magnitude relative to the current measurement. In a few years, 

—)■ 7t~^u17 will be one of the the most sensitive havor observable to Q-class 
models. Perhaps more importantly, if NA62 were to measure a deviation from the 
standard model prediction, this model would provide a signihcant motivation to 
invest in an experiment like ORKA [66], that would be able to hone in on the 
parameter space. 
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Observable 

Projected Accuracy 

ArriK 

10%,,. [67] 

Ams, 

10%„. [61] 

AniB, 

5%„. [61] 

Amo 

None 

Br{K^ —)■ 

10%e.p. [61] 

Br{B ^ X,7) 

7%e.p. [63] 

Br{B ^ X,7) 

24%,.,. [63] 

Br{Bs -)■ /iV) 

15%exp. [61] 

Br{Bd /iV) 

35%e.p. [61] 


ATi/a-j 

Figure 9: Future projections for flavor constraints on EGMSB in the ~ 3-5 year range. 
In the hgure, the errors have been updated to the fractional values in the accompanying 
table, and the central values have set equal between theory and experiment. Depicted is 
the same Q-class model from hgure 8. The dashed green line, which is unchanged from 
hgure 8, indicates where the EGMSB contribution to ArriD is equal to the measured 
value. Observables not shown are not constraining. 

Due to these potentially signihcant theoretical and experimental improvements, much 
of the parameter space in Q-class models could be probed in just a few years. The pro¬ 
jected sensitivities are shown in hgure 9. Although Q-class models have an exciting 
future in havor, t/-class models remain completely unconstrained. Short of lattice pre¬ 
dictions for Amo, no observables in the current program are sensitive to these these 
models. Gharm factories and precision top studies could someday explore this space in, 
for instance, c —)■ uy or f —)■ c/uj. However, for the moment, [/-class models are resilient 
against havor constraints. 

4.3 Future Directions 

There are many avenues for future investigations. Here we list a few. 

• As we alluded to earlier, the ansatz is a general paradigm, and it provides a 
novel and realistic solution to the SUSY havor problem. This texture, its possible 
origins, and consequences are worthy of further study [51]. 

• Another interesting question is that of GP violation in these models. The goal of 
































this work was to focus on the SUSY flavor problem, and so the EGMSB couplings 
were intentionally assumed to be real. Allowing for arg k* 7 ^ 0 would give nontrivial 
contributions to CP-violating observables, some of which (such as ex and the 
neutron EDM) are typically extremely constraining. An interesting question is 
to what extent the texture protects EGMSB models from the SUSY GP 
problem before there is any conflict with data. We plan to study this in detail in 
an upcoming work [35]. 

• The possibility of heavily mixed squarks allows for very interesting collider signa¬ 
tures [68-71]. These EGMSB models provide a flavor-safe proof-of-concept moti¬ 
vation for experimental searches at ATLAS and GMS. 

• While we made some effort to ensure that the flavor-violating A-terms do not 
destabilize the vacuum (see the discussion in appendix A), it would be interesting 
to study in more detail the vacuum stability of these flavor-violating EGMSB 
models, along the lines of [72]. 
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A Details of the Deformation and Parameter Scan 

In this appendix, we detail our parameter scan. In particular, we specify the deformation 
about the best points in [20] that we use. Our type I models are characterized by the 
following parameter space: 


(ki, K2, K3, A/^,A) (A.l) 

In [20], the models were studied at the Ki = K 2 = 0 point. For each % and A/M, A 
(which sets the overall scale of the superpartner spectrum) was increased until mh = 
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125 GeV was achieved. Then the hne-tuning of the point was estimated and the re¬ 
gion of least hne-tuning in (k 3 ,A/M) space was identihed. Our aim is to investigate 
these regions in the presence of nontrivial havor violation, parametrized by ki and 
K 2 - a full optimization of the hne-tuning would involve a hve-dimensional scan in the 
(ki, K 2 , M) parameter space. Such an endeavor is not computationally feasible 

and moreover is unnecessary. The main question we intend to explore is how qual¬ 
itatively dangerous the havor-physics contributions are in these EGMSB models. In 
principle, this question can be answered using a simpler scan where we essentially hx 
(^ 3 , A/M, M) as in [20] and then vary k\, K 2 - However, there are a number of subtleties 
to take into account when doing so, most of which stem from the fact that the points 
from [ 20 ] tended to have light stops as a result of a cancellation (this is unsurprising due 
to the tuning involved). 

• We must be careful to choose {k 3 ,A/M) such that the point at ki = ^2 = 0 is 
unconstrained by Run I searches. In practice, to get a model that is not constrained 
by the LHG and has a good tuning value, we will lower A away from the least-tuned 
point identihed in [20]. This increases the hne-tuning required in the models only 
by about 10%. The parameters we choose (at the origin) are K 3 = 0.858, N = 6, 
^ = 0.347, and M = 312 TeV for 1.9, and K 3 = 0.908, N = 3, ^ = 0.290, and 
M = 350 TeV for 1.13. The spectra are shown in hgure 1. 

• Despite our focus on GP conserving observables in this work, there could, in prin¬ 
ciple, be signihcant constraints from eK introduced solely through the GP violating 
phase of the GKM.^° In order to avoid this constraint, in our Q-class models, we 
choose Ki, K 2 and K 3 to align with the low-energy down, strange and bottom quark, 
respectively. Thus in our choice of interaction basis, the down Yukawa is diago¬ 
nal, but the up Yukawa is multiplied by Vckm- We ignore the small diherences 
with our previous work where the alignment was with the top quark. The [/-class 
models retain alignment with the up-type quarks. 

• Turning on ki and K 2 while holding hxed the other parameters can signihcantly 
modify the spectrum, as shown in ( 2 . 6 ). At the least hnely-tuned points, where a 
cancellation results in a stop lighter than the other squarks, this can either lead 
to stop tachyons, or it can lead to the stops being so heavy that the hypercharge 
tadpole contribution to the RG running quickly drives the sleptons tachyonic. To 
avoid these undesirable features, as we turn on 2 , we fix A, but vary M, so that 

thank W. Altmannshofer for bringing this to our attention. 
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the the one-loop contribution adapts in such a way that the lightest squark eigen¬ 
value of the LL block in Q-class or RR block in [/-class models is held hxed. This 
deformation enables us to prevent the interesting flavor-violation from vanishing 
and to maintain the squark masses at sensible values. 


• Additionally, as can be seen from (2.6), if we were to £x K 3 and turn on ki^ 2 , then 
the “net” A-term would increase. The “net” A-term in the type I models is aligned 
with the lightest eigenvalue direction of the rank 1 block, i.e., in Q-class models. 


Cd- 




HuQstR- 


(A.2) 


The A-terms in all orthogonal squark directions vanish in the third-generation 
dominant limit. Increasing A/ms is potentially dangerous for vacuum stability 
[72]. In order to avoid these issues, we require that this “net” A-term remains 
constant. To achieve this we £x k^k everywhere in the parameter space so that as 
we increase ni and K 2 , decreases to compensate. 


• Finally, this deformation of ^ 1 ^ 2,3 will also affect the Higgs mass. The one-loop 
corrections to the Higgs mass in the presence of general flavor violation have been 
computed in [73, 74], and our preliminary studies of these suggest that may drift 
down by several GeV as we move out in ki 2 . However, the two-loop corrections 
are not yet known, and if the usual non-flavor-violating MSSM is any guide, these 
are likely to be important for an accurate determination of the Higgs mass. While 
it would be interesting to study this further, it is beyond the scope of the work. At 
the very least, one could imagine increasing the overall scale of the superpartners 
in order to compensate for any decrease in the Higgs mass. This would only serve 
to further alleviate the flavor constraints, so the qualitative value of XFV to these 
EGMSB models is unaffected. 


Across a grid in {ki, K 2 ) constructed via this deformation of the benchmark point, we 
generate the soft spectrum at the messenger scale. All of the couplings and soft masses 
are then evolved down from the messenger scale to the SUSY scale using the fully general 
3x3 MSSM /9-functions.^^ At the SUSY scale, we apply the BMPZ QGD threshold 
corrections to the squarks and gluinos [75]. Finally, we use our new Mathematica package 
FormFlavor to compute the Wilson coefficients, RG evolve these coefficients down to 
the scale of interest (e.g., m^) for each flavor observable, and compute the contributions 

^^For simplicity, we do not run back and forth between the IR and UV to better specify the scale. 
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Ki/k 3 = ^2/^3 


Figure 10: The spectra of squark masses (in GeV) for model 1.9 as we move in ki = K 2 
Aside from the three lightest eigenvalues, ms^up = ^s,down, and mjui, the other nine 
squark eigenvalues all fall in the range of the blue band near 3.5 TeV. The width of the 
band is largely due to the SU(2) GMSB contributions to the left-handed squarks. 

to each flavor observable there. (The details of FormFlavor are briefly explained in 
appendix B and will be further fleshed out when the package is made publicly available 
[39].) 

B Details of FormFlavor 

FormFlavor is a general, modular flavor package written in Mathematica that takes 
general MSSM spectra as inputs and calculates a variety of flavor observables in situ 
starting from the Feynman rules. The package is flexible - new flavor observables can be 
straightforwardly added to it and, although currently implemented only for the MSSM, it 
can be extended to other models readily. Also, since processes are calculated in a uniform 
manner from hrst principles (rather than hardcoding formulas from the literature), the 
reliability of the code is greatly enhanced. 

In FormFlavor, the process under consideration is specihed in FeynArts [36] and all 
the relevant topologies for the various sparticle mediators are generated there. These 
topologies are then converted into their respective amplitudes via FormCalc [37]. The 
FormCalc output is analytically transformed into a particular Wilson operator basis for 
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each process of interest, and the Wilson coefficients are extracted. Having performed 
this once for an observable the analytic expressions may be written to a process hie and 
used in subsequent runs. The Wilson coefficients may then be piped into functions that 
calculate the havor observables. 

The main components of FormFlavor are as follows: 

• Automated generation of one-loop diagrams for havor processes (using FeynArts) 

• Automated calculation of Wilson coefficients from Feynman diagrams (using FormCalc) 

• Library of analytic one-loop integral functions 

• Routines for general MSSM spectrum input (SLHA2 [76] compatible) 

• Routines for converting Wilson coefficients to havor observables. (This step con¬ 
tains hardcoded formulas from the literature and is not MSSM specihc.) The 
following processes are currently implemented: 

CP-conserving observables 

* Meson Mixings [52, 77]: AM^ [78], AM^ [79], AMb^ and AMb^ [80, 40] 

* Leptonic decays : —)■ [81, 41, 82], Bg —)■ and —)■ 

[55-57, 46] 

* Radiative Decays : b ^ s'j [83, 42], 6 —)■ d'j [84, 44] 

CP-violating observables 

* Meson Mixings : ex 

* Leptonic decays : —)■ 

* Electric dipole moments : neutron-EDM 

• MSSM RGE - Full three-family MSSM renormalization group evolution, including 
arbitrary CP-phases, implemented with the aid of SARAH [85] 

This powerful havor package will be released to the public and described in more detail 
in an upcoming work [39]. 
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C Soft Terms for Type I Squark Models 


In this appendix, we will specialize the fully general EGMSB formulas from [20] to the 
case of the type I squark models with general flavor-violating couplings Ki. These were 
quoted in (2.6) for the Q-class models, 

N 

wEGMSB = l^iQi ^ (C-1) 

A=1 

and we again focus on this case here. To translate to [/-class models, one must take 
U y-)- Q, ?/„—)■ j/1, As —)■ At and K -A- in all expressions. 

Before we begin, we must address a convention difference between this paper and 
that of [20] concerning the treatment of the right-handed quark chiral superhelds [/, D. 
In this paper, the Yukawas, squark mass matrices and A-terms are dehned as: 


hh Z) H^Uiyy^ ijQj HdDiyfj^ ijQj 

C D (C.2) 

+ {HuU* Au^ijQj + HdD*A^-jQj + h.c.^ 


Thus, the Yukawas, A-terms, mfj and are all transposed relative to those of [20]. 
Note that in these conventions, U*, D* are the scalar components of the chiral super- 
multipliets U, D. 

The starting point of our derivation of (2.6) is eq. (2.19) in [20], which we repeat 
here for convenience: 


Aah = 


KsC^bBC^ 


^'^ab — 25g7j-4 ^'b^^ aBC^bCEl^cBol^cDE + X^^fjXcBC^cDE^ 


“0 ■-aBC'^cBcXcDEXbDE 


(C.3) 


“ 2 ^a “c yacdUbdeAcBC^eBC ~ “a 9r^aBC^bBC I A 


Here a,b, ... run over all MSSM helds (including all gauge and flavor degrees of freedom). 
Meanwhile A, B, ... run over all messenger helds similarly. In this appendix, we will 
use i, j, ... to denote MSSM havor indices. 

From (C.3), we immediately obtain the bilinear A-terms after summing over the N 
messenger multiplets and substituting Xasc k* and df*" -A- dg oc N 


A 


Q,ij 


1 

lOvr^ 


dQK*KjA 


(C.4) 
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This becomes the trilinear A-term used in the paper via: 


yu,ikAQ^kjHuU*Qj = Au^ijHJJ*Qj 


(C.5) 


Next let’s consider the EGMSB contributions to the soft mass-squareds, starting 
with ttIq. From (C.3), we obtain: 


_ _ I _ 

- 2567r4 




Q 


2567r4 


dQd^K*KjKkKl + dgnluluki^j - 2dQCrg^tA^Kj A 


{d^ + dQ)K^ - 2Crgl 


(C.6) 


In the hrst line, we have introduced djP —)■ ^d^, and we have used the fact that Q is the 
only MSSM held coupling to the messengers to set d^^d^^ to zero. Additionally, as Q 
couples directly to the messengers, there is a one-loop term suppressed by [13], 


5 (m 




Qm I 


M) 2567r"^ 


(C.7) 


where h{x) is a loop function given by 

~ 2) ln(l — x) — {x + 2) ln(l -|- = 1 -|- -|- O (x^) . (C-8) 


The EGMSB contribution to is much simpler. Here only the third term of 


ab 


in (G.3) contributes: 


Smir = 


3d, 


Q 


3dQ 


2567r4 

where we have used d^^ = 3. 


Tr [yuKyi] A^ 


(G.9) 


Lastly, the EGMSB contribution to also comes from just the third term of 




d?r^d, 


u “Q 


y*u,ikyjif^kK*iA‘^ = 


d?r^d, 


u “Q 


{yuKyi) A^ 


(G.IO) 


2567r4 ^ 2567r4 

where d^^ = 2. Taking into account the need to transpose to translate between the 
conventions of [20] and those of this paper, we obtain the correct result quoted in (2.6). 


D Formulas for flavor observables 

In this appendix, we collect formulas from the literature for the various havor observ¬ 
ables considered in this work. Along the way, we will streamline the different notations 
scattered throughout the literature into a uniform convention. 
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The uniform operator basis we will use was introduced in the text; we repeat it here 
for convenience. For dimension 5, we have: 


0^{hj2) = ef,a'^^PMf2F,. 
0^{hj2)=gfiCJ^'^PMf2G,. 


(D.l) 

(D.2) 


For dimension 6 , we have: 


/2. A. A) = CfiPMh)QzPNU) 

A. A) = (fii-PMf2)Qii^PNf>) 
of'"(A. A, A. A) = (7iff'"'ft<A)(7sff;.-^’A'A) 


(D.3) 

(D.4) 

(D.5) 


where M, N = L,R. Pr = ^ (1 + 75 ), | (1 — 75 ) are projection operators, = 

^[ 7 ^, 7 '^] and if /j carry color indices, they are contracted within a bilinear factor. 

The general effective Hamiltonian is then: 



(D. 6 ) 


where the sums runs over a complete basis of independent operators. 

D.l Meson Mixing Anix 


Meson 

mx(GeV) /x(GeV) 

Rx 

nLL 

By 



tjLL 

P 5 

tjLL 

hjrp 

^^y,WU=-,(GeV-2) 

ArriK 

0.4976 

0.160 

24.3 

0.56 

0.85 

1.08 

0.62 

0.43 

- 

Amo 

1.8645 

0.209 

3.20 

0.76 

0.97 

0.95 

0.64 

0.39 

- 

Ams, 

5.2796 

0.191 

1.65 

0.84 

1.47 

0.95 

0.72 

0.61 

(2.34- 2.20i)xl0“^2 

Ams, 

5.3668 

0.228 

1.65 

0.88 

1.57 

0.93 

0.73 

0.62 

(6.96 -0.26f)x 10-^1 


Table 2: Properties of the meson mixing used in this work, fx and Pyx come from 
the FLAG review [77]. The other four non-perturbative P-parameters are taken from 
several sources: for ArriK from [78] (at /i = 2 GeV), Amn from [79] (at /x = 3 GeV, 
rescaled for a common Rd, and converted to our basis using = |P 2 — fPs), and 
for Bs and Bd from [80] (at ^ = rrib = 4.2 GeV and converted to our basis). Cy^j^\f^=mf, 
are the FormFlavor values in the GKM basis of the PDG [47]. 

This corresponds to a AF = 2 effective Hamiltonian, with fi = fs = qi, and 
f2 = fi = ^ 2 , with (gi, ^ 2 ) = (s, d), (c, u), ( 6 , d), ( 6 , s) for K, D, Bd and P^ respectively. 
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The quantities relevant for mixing are derived from the effective Hamiltonian as 


{X\H,„\X) = Mx,,2 - ^Tx,u 
{X\H,f!\X) = Mi, 12 - Ti,i 2 


(D.7) 


Mx,i 2 and Tx,i 2 are, respectively, the dispersive and absorptive parts of {X\Heff\X). In 
terms of these quantities, the mass splitting is given by: 


Xmx = 2 Re 


Mx,i2 — 


( M* — IF* 
^^■^X,12 2 X,12 

Mx,i2 — |rx,i 2 


(D.8) 


For X = Bd and Bg where rx,i 2 ^ Mx,i 2 , this is well approximated by 

Amx^2\Mx,i2\-2\{X\H,ff\X)\ 


(D.9) 


For X = K and D, where experimentally one finds that Mx,i 2 /^x,i 2 is approximately 
real, and Mx,i 2 and Fx,i 2 are both predicted to be approximately real in the standard 
CKM convention (where the CPV phase is primarily in Vtd and Vub), one has to a good 
approximation 

Amx ~2ReMx,i2 ^ 2Re {X\Heff\X) (D.IO) 

Finally, the short-distance part of the matrix elements in (D.9) and (D.IO) (which is all 
that is relevant for the Bg systems and for new physics) are given by 


{X\H,ff\X) = 


mxfx 

24 


- Rx 


ABb^Ch^ - 6Bb^C 


LR^LR 


'V 


s 


+ + 12 R|;^ 


(D.ll) 


where. 


Rx — 


mx 


rug^ m, 


92 


(D.12) 


and the R-parameters are non-perturbative corrections that have been computed on the 
lattice. In the Standard Model, only ^ 0. These values and the other parameters 
relevant for meson mixing are shown in Table D.l. 
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D.2 




This corresponds to a dimension 6 effective Hamiltonian with /i = s, /2 = d, /s = /4 = 
z/£. The branching ratios for —)■ are given by [81], 


BR(A'± ^ = ^1,“ ^ + 

i=e,ii,T 


(D,13) 


where v = 246 GeV and 


c+ 


3rj^+ 

2|K.| 




1.35. 


(D.14) 


Here, the branching ratio has been inclnded to remove dependence on the hadronic 
matrix element, and r^+ = 0.901 contains isospin violating qnark mass effects and 
electroweak corrections compnted in [86]. 

In the standard model, the contribntions from top loops are have no sensitivity to 
different generations, bnt the charm loop contribntions do. In particular, the charm 
contributions are the same for e and /r, but differ for r. Thus, the SM contribution can 
be expressed as, 

CfsM/= ^ (Wf + A,X,) (D.15) 

where = V*Vid, Xt = 1.469 [41], = 1.055 x lO-^, and X^ = 7.01 x 10"^ [81]. 

In the absence of lepton flavor-violating new physics effects, the contributions from 
new physics are the same across lepton generations, i.e., CyJ^pf^ = CyJjpp = CyJjp. 
We can then give a simpler form to equation (D.13) 


BR(X± ^ TT^i/P) = c+v^ 


(D.16) 


where 


r'XY I r'XY 
'^V,NP “T 


with 


^RL _ 

'^V,SM — 

p^LL _ 
'~'V,SM — 


Pc = 


0 


^ (AePc + AiXt) = (-12.2 + 3.6i) x 10“^^ GeV^ 

+ G; A ~ 9.37 X 10-. 


(D.17) 


Here, Pc is the charm contribution averaged over the different neutrino flavors. Using 
Pc simplihes the expression, but reduces the standard model charm contribution by 
about 3%. This represents only a 0.3% decrease in the overall SM contribution, which is 
completely negligible when compared to the theoretical uncertainty. Most importantly, 
the interference effects with new physics are properly captured under this simplification. 


38 





D.3 


These observables correspond to dimension 5 effective Hamiltonians with fi = b and 
/2 = s or d. The branching ratio is given by, 


where [87] 




= 

$ = 


(Svr^)^ 

Vub " 
Vcb 


6 BR(6 —)■ Xcez/)Bxp 
TT $\Vcbf ml 

BR(6 -»■ XcCu) ^ „ c-o 
BR(6 ^ X^eu) 


Q^EM 


3.3 X 10 ^ 


(D.18) 


(D.19) 


where the branching ratio to charm decays is nsed to remove sensitivity to the hadronic 
matrix element, and the <h factor is introdnced to acconnt for the nontrivial phase space 
factor in the compared branching ratio (due mostly to the charm quark mass). 

Only contains a standard model contribution, which is. 


~ (1.3 + 0.03i) X 10-* GeV (D.20) 

where Xsm = —0.3736. 

The observable 6 —)■ dy is defined completely analogously. The only difference is that 
the the standard model contribution differs. In particular, 

~ + !■«) X 10-’GeV-‘ (D.21) 

with the same Xsm = —0.3736 as above. 


D.4 Bg ^ /i+/i and 

This corresponds to a dimension 6 effective Hamiltonian with fi = b, /2 
/s = /4 = /^- The branching ratio for Bg d —)■ /i+p,“ is [55, 56]: 

BR(R, ^ /iV) = X, I (^1 - 

where 


s, d and 


(D.22) 





1287rmsTB^ 



X., 


5.36 X 10 ^ and X^ = 3.97 x 10 ^ 


(D.23) 


39 













(D.24) 


= 


m 


B, 


TUb + m,, 


(/^LL I ^LR ^RR /^RL\ 

A'-^S +*^5 *^5 '-^5 


m 


B, 


(-C|^ + 


RL\ 


= 

nib + Bfii 

Ff = 2mB,m^(C'“ - 


and = 4.33 x 10 GeV and = 4.49 x 10 GeV. The three-loop 

Standard Model contribntion to Bi —)■ can be expressed as [57]: 

,2 4 ^ksM ~ ^p!sm ~ 0 (D.25) 


F%m = O-^SOS X “““a ■ >^“1 - 


TT 

SO that F^gj^ = (1.5 — 0.6^) x 10“® and F^^gj^ = (—7.9 — O.li) x 10“® 


E Rank 1 xFV Loop Functions 

We compile some of the loop fnnctions for our rank 1 X^V approximation below. 
Starting with meson-mixing, we have, 


i^AA/.box, 

fa Fg,Xq) = 


(E.l) 


2Xq loga:g(llXg -|- GXgX^ — 2x^Xq — ISXgXq — 2x^) 

{Xq-l){Xq-XgY 

2XqXg log Xg[Xq - l)'^ (2XqXg + 13Xq - 17Xg + 2Xg) 

{Xg-lf{Xq-Xgf 

Xq (l9XgXg + lla:^ -f- 30:^0:^ — 74xqXg -7 llXg -f- 8Xg + 3Xg -f 19xg) 

{Xg-lY{Xq-Xgy 

Here, Xq = m^/nil and Xg = my ml] where mo and ms are the heaviest and lightest 
squark mass eigenvalues. 

For —)■ TT^z/I/, there are two contributions. 


(E.2) 


fy-^-^’^{xe,X2,x,) = fy 

In practice, 0 : 2 ) is numerically only 0(10%) of 0 : 2 ), so, out 

X X 

of simplicity, we only present the box contribution for interpreting our results, 

X 2 (x^ - Xe) logXs_X 2 _ 

2 (x2 - 1) fx2 - Xf) 

(E.3) 


K^7riyiy,hox 


4± 


(X(>,X2) = 6 


(X 2 - 1 ) 2 (x 2 - 


(X 2 - 1 ) (X 2 - Xi) 
X2X1 logX£ 

(X2 -xf)‘^{xi - 1 ) 
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Here xi = X 2 = M^/m^g and x^ = /m^g. We have also dropped the neutralino 

boxes and Z-penguins that contribute at the O (10%) level or below. 

For b ^ q'j 


^6^s/d7,peng / 


, 3 : 2 ) 


n 


5x2 Xfj_ — 7x^ — 7x2 + 9 2 (2x2 — 3) log X 2 

(X2-I)'(x^-1)' ^ (a;2-l)=*(a;2-a;^) 

2(2x^ - 3 ) logx^ 

(x 2 - x^){xf, - 1)3 


(E.4) 


Here x^ = ^‘^/rrig, X 2 = M^/rrig and ms is the lightest up squark mass eigenvalue. 
For Bq —)■ 


pBq—^ll+fl ,h-peng/ 

Jg {Xq, 


Xg) =y/XgXq 


(Xg - l)XglOgXg 1 

(Xg - 1)^ (Xg - Xg) {Xg - l) 


Xq log Xg 

(Xg - 1) (Xg - Xg) 


(E.5) 


nBq—>-g+ll ,h-peng/ 

Jx \Xq^ 


X 


gi 


X 2 ) =y/^2 



X 2 (Xg - 1 ) log X 2 

(X2 - 1 ) (X2 - Xg) (X2 - Xg) 

Xg{Xq - 1 ) lOgXg _ 

- X 2 ) (Xg - 1 ) (Xg - Xg) (Xg 


Xq log Xg 

- Xq) (X2 - Xg) 


(E.6) 


Here Xg = Xg = jj?X 2 = M^/m^ and and Xg = m^g/m^] ms and mg in this 

context are the lightest and heaviest squark mass eigenvalues. 
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